Exposing the non-collectivity in elliptic flow 
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We show that backward-forward elliptic anisotropy correlation provides an experimentally ac- 
cessible observable which distinguishes between collective and non-collective contributions to the 
observed elliptic anisotropy V2 in relativistic heavy ion collisions. The measurement of this observ- 
able will reveal the momentum scale at which collective expansion seizes and where the elliptic 
anisotropy is dominated by (semi)-hard processes. 
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Deconfined QCD matter at high energy density, the so- 
called Quark-Gluon Plasma (QGP), was a phase during 
the evolution of the early universe and it is now created 
and explored experimentally in relativistic heavy ion col- 
lisions. One major discovery by the experimental pro- 
gram [l[ of the Relativistic Heavy Ion Collider (RHIC) 
is the large elliptic anisotropy Q, of observed parti- 
cles' transverse momenta, which is consistent with predic- 
tions from ideal hydrodynamic expansion 0] . This led to 
the conjecture that the matter created in these collisions 
exhibits "perfect fluidity" , i.e. minimal shear viscosity. 
First studies 0, Q indicate that the shear viscosity 
must be very small especially near T c , and likely smaller 
than any known condensed matter substances and rather 
close to the conjectured universal lower bound 2 > ^ 
based on calculations utilizing the gauge/string duality 

0- 

The transverse momentum (p t ) dependence of the el- 
liptic anisotropy, u 2 , as depicted in Fig. Q] shows a rise 
at low pt towards a maximum at p t ~ 3 GeV, and a 
constant value for large transverse momenta. At present 
the rise at low transverse momentum (p t < 1.5 GeV) is 
thought to be due to collective hydrodynamic expansion, 
which translates the initial spatial anisotropy into a pt 
anisotropy over a wide region in pseudo-rapidity. 

At very high p t , on the other hand, is believed 
to result from the different attenuation of hard partons 
(jet quenching) in the asymmetrically distributed mat- 
ter [3, 1^, 11 1 . In this case the elliptic anisotropy, V2, is 
non-collective in the sense that it is due to jet-like pro- 
cesses which are rather local in rapidity [12|. The tran- 
sition from collective (hydrodynamic) to non-collective 
(jet-like) anisotropy is expected to take place around 
Pt — 2 — 4 Gev, but the details are not well understood. 
In addition, one would expect that the attenuation of the 
jets should result in local, non-collective anisotropy also 
at lower p t , since the debris from the jet-quenching pro- 
cess needs to go somewhere. If there is indeed a sizeable 
non-collective contribution to vi at low p t , comparisons 
of (viscous) hydrodynamics with the data may lead to 
wrong conclusions about the viscosity of the produced 
matter in these collisions. Therefore, it would be desir- 



able to have a direct measurement not only of the pt scale 
at which the transition from collective to non-collective 
vi occurs, but also of the contribution of non-collective 
effects in the low p t region, p t < 1.5 GeV. 

It is the purpose of this work to propose an observable 
which will distinguish between the collective and the non- 
collective contributions to the elliptic anisotropy. In this 
paper, for illustration purposes, we will identify the col- 
lective component with hydrodynamic flow and the non- 
collective one with attenuated jets. We note that there 
may be other mechanisms at work which generate sim- 
ilar collective and/or non-collective contributions to the 
elliptic anisotropy, and the proposed observable is not 
sensitive to any specific mechanism. 

The key observation is the following: consider two ra- 
pidity bins, one forward, one backward, with a suitable 
separation (gap) in between. In a given event, the elliptic 
anisotropy of the collective component in both rapidity 
bins is highly correlated. The non-collective component, 
on the other hand, generates an elliptic anisotropy in ei- 
ther the forward bin or the backward bin, but never in 
both, resulting in very small backward-forward correla- 
tions. To illustrate this point, consider hydrodynamic 
expansion as an example for the collective component. 
In each event the elliptic anisotropy is aligned in az- 
imuth and of comparable magnitude in the forward and 
backward bins. Contrast this with the elliptic anisotropy 
due to quenched jets, as an example for a non-collective 
component. At RHIC energies one rarely has more than 
one hard process per event. In addition, the anisotropy 
due to attenuated jets will be very local in rapidity and, 
therefore, will cither contribute to the forward or to the 
backward bin but not to both. Consequently, these pro- 
cesses do not give rise to long-range forward-backward 
correlation of the elliptic anisotropy. At higher energies, 
such as the Large Hadron Collider (LHC) with many hard 
processes per event one has to choose a suitable rapid- 
ity gap, which is larger than typical jet cone extension 
but smaller than average jet-jet separation. After these 
general considerations, we will proceed to define the pro- 
posed observable and demonstrate its sensitivity to the 
degree of collectivity in contributions to V2- 
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FIG. 1: Parametrization of V2 data for AuAu 200 GeV, (00-20% centrality). The red diamonds are for negative charged 
hadrons [ItJ and the green boxes are for no from PHENIX Run7 Preliminary fl3j fill]. The blue/orange/magenta solid lines 
in (a)/(b)/(c) represent the combined parametrization, Eq.@, respectively. In all three panels, the black solid lines show the 
collective-flow-only contribution while the colored long-dashed lines (extending to high p t ) show the jet-only contribution in 
Eq. (|13[l , The dash-dotted lines show the weighted contributions (1 — g)- < > (black) and g- < v% > (colored), respectively. 



Let us start by recalling the definition of v 2 : 
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is the distribution of the event- yield and the jet-related yield, respectively: 



In the above x — , — ,- 

\ ptdptdtp 

averaged ^-differential yield over the azimuthal angle 
<f> which is defined with respect to the reaction plane. 
The numerator, denoted as V2, may be referred to as 
the total elliptic anisotropy. Here we assume that the 
reaction plane is determined with high accuracy as it 
is achieved by current RHIC experiments. The particle 
yield at RHIC can be attributed to two main sources: the 
bulk matter which dominates low pt regime and exhibits 
collective flow, and the (partially suppressed) hard jets 
dominant at high p t . Both sources contribute to the total 
yield d , N ,> and thus to the measured v 2 . 

J ptdptdcfi z 

The proposed new observable Cfb[pt] is the correla- 
tion of the total V2(pt) between forward (F) and backward 
(B) rapidity bins, specifically 



[±y m i n , ±y max ] should be chosen such that the gap 
{—ymin, ymin) is sufficiently wide to prevent a jet from 
contributing to both simultaneously while still keeping 
both bins within the "flat plateau" near mid-rapidity. 

The total yield, = + can be decom- 

posed into two components [Ioj|. the hydrodynamic flow 



Cfb [pt] = 



<V 2 F -V2 B > 

< Vo F > ■ < vf > 



(2) 



where the total elliptic anisotropy in the forward, V 2 F , 
and backward, , are defined as 

/ \ r 2rr /hn f / b \ 

(v 2 F/B ( P t)) = j a d^cosmi-^y (3) 
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(V 2 F -V 2 B ){ Pt ) = {J d^cos(20)^ 
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Here we denote by = ^ ^ ■ Pt ■ Ap t the 

differential particle yield for a given p t interval in the 
F/B bins, respectively. The forward/backward bins 
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(5) 



F+B 



(In the following we drop the a F + B" subscript.) 

To quantify the main idea, we express the yields in the 
F j B bins in each event as: 
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In the above we have introduced three random variables 
to schematically describe the fluctuations from event to 
event. f] F ^ B represent (independent) random deviations 
from the average hydro flow yield in F/B bins, and satisfy 
< T] F I B >= 0, < r] F ■ i] B >= 0. The variable £ assumes 
values of either or 1 in each event with equal probabil- 
ity, i.e. < i >=< 1 - f >= 1/2 while f (1 -0 = 0. The 
physical idea is that in each event there is at most one 
high p t hadron cluster contributing to the final observed 
anisotropy V 2 [13] either in the forward or the backward 
rapidity bin. Calculations of dijet production at RHIC 
energy show a rapid decrease of cross section with in- 
creasing rapidity separation [lij. Also with a hard trigger 
in one bin, the back jet is strongly degraded in a heavy 
ion collision, inducing negligible contribution to v 2 . In 
contrast, the hydrodynamic flow (ideal or viscous) con- 
tribution to the anisotropy is about equally split between 
F/B. The so-called anti-flow [l5j gives rise to correlations 
between rapidity y and in-plane p x . Its magnitude, i>i, 
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at RHIC energies is at most a few percent 0, [l5[ and in- 
duces a correction to t>2 of the order Sv2 = v\ < 1%. 
Furthermore it contributes equally to vi in both the F/B 
bins due to its (anti-)symmetry with respect to y = 0. 
Next we introduce an interpolation function: 
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with g(pt) and 1 — g(pt) giving the relative weight of the 
jet (non-collective) and hydro flow (collective) contribu- 
tion to the total, F + B, yield, respectively. Given this 
function we express the observed < v% > in the F + B 
bins via ([!]) in terms of the flow and jet contributions: 



< v 2 (pt) >= 



r#cos(20)[<^f) + <^f)] 



= [1 - g{p t )\ < wf > +g(p t )- < v{ 



> 



(9) 



with < v^ jJ > defined as in (TT]) with the corresponding 
hydro-flow-only/jet-only yields. Combining (13I4I6I7I8I9[) 
and using the formalism of [l6| . we find for Eq.([2]): 



(1 - gf < vl > 2 +2g (1 - g) < >< 
Cfb 77 —7^772 



> 



[(1 - 9) < «f > +g < vi>y 



(10) 



A distinct feature is that as g — > (hydro dominance) 
Cfb — * 1 while Cfb — > in the other limit, j -> 1 (jet 
dominance). Thus, the correlation Cfb [pt] distinguishes 
between collective and non-collective contributions and 
exposes the transition from the flow to the jet regime. 

In order to demonstrate the sensitivity of Cfb we 
study three different scenarios within a simple model. 
All the scenarios employ the same blast-wave model for 
the collective flow contribution, but they differ in the pt- 
dependence of the jet-only elliptic anisotropy, (v£ (pt)), 
which is not known for small p t < 3 GeV. The rela- 
tive strength, g(pt), of flow and jet contributions is ob- 
tained by fitting the V2 data of the PHENIX collabora- 
tion [H, [l7| El for AuAu - 20% centrality class at 
V^s = 200 GeV (see Fig[T]) . Scenario (a) assumes a con- 
stant (v£ (pt)) (see FigUta)) . Scenarios (b) (FigQTb)) 
and (c) (Figfljc)) assume a constant (v$ (j>t)) for large 
p t which drops to half its value at pt = for scenario 
(b) and to zero for scenario (c). For both scenarios the 
momentum scale at which (v£ (pt)) changes towards the 
values at low p t is pt = 2 GeV. 

The blast-wave model for the collective hydro flow con- 
tribution we adopt from ref . [191 ] . By parameterizing the 
flow velocity field at freeze-out as in [19( we have, 
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The elliptic anisotropy is intrinsically built in the flow 
field u^(x^): in (T,r/ S ,r, <fi s ) coordinates u r = u [l + 
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FIG. 2: The interpolation g(pt), Eq. (114[l . for scenarios (a) 
(solid line), (b) (long-dashed line), and (c) (short-dashed line). 



u 2 cos(2(j) s )} 0(R o - r),u T = y/l + (u r ) 2 , = = 0. 
The blast-wave parameters are chosen as = 140 MeV, 
T = 170 MeV, R = 10 fm, r G = 7 fin, u a = 0.7, and 
U2 = 0.06 which approximately reproduces the yield at 
low pt [l3J] . Given this model we calculate the flow-only 
< V2 > via (TTJ using the yield in (fTT|) (black solid lines 
in FigQJa,b,c)). We also show the weighted flow contri- 
bution (1 — g) ■ < v% > as the black dash-dot lines. 

In order to model the anisotropic jet attenuation we 
write the jet-related yield as 

(^~) CX [l + 2<«f)(pt)-cos(2#] (12) 

where < v£ > (pt) parameterizes the resulting az- 
imuthal anisotropy. While the experimental data indi- 
cate < i>2 > (pt) to be sizeable (7%) and rather constant 
for p t > 6 GeV [l8| , it is not clear how v£ behaves at 
lower p t . To explore this we study the three scenarios 
described above which we parameterize as follows: 



< vi >= 0.07 [a + (1 - a) tanh(p t - 2 GeV)] (13) 

The three scenarios (a), (b), and (c) correspond to a = 1, 
a = |, and a — i, with the respective < v?f > (pt) 
plotted in FigQ] as colored long-dashed lines. 

Finally, we parameterize the interpolation ([5]) as: 



g(pt) = [1 + tanh[(pt - P c )/Pw}]/2 



(14) 



The parameters Pc , Pw ( m units GeV) are obtained by 
X 2 fitting of Eq. ([9]) to the data for V2 ■ The best choices 
for the three scenarios are: (a) Pc = 1.4, Pw = 2; (b) 
P c = 2.4, P w = 1.8; (c) P c = 2.4, P w = 1.8. In all 
cases a reasonable model description of the v 2 data is 
established over the whole pt range, as can be seen in 
FigfT] The resulting interpolation g(pt) are plotted in 
Figf2] At low to intermediate p t they differ consider- 
ably in the magnitude of the non-collective component. 
Of course the plotted curves are only examples for what 
may happen. To quantify g(pt) more experimental infor- 
mation will be required. We note, that an admixture of 
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FIG. 3: The correlation Cfs(p*) for scenarios (a) (solid line), 
(b) (long-dashed line), and (c) (short-dashed line). 



transition pattern would be of interest, as the flow and jet 
components scale differently with centrality. In addition, 
these measurements will help to understand the origin of 
the decreasing vi at intermediate pt and constrain viscous 
corrections to the hydrodynamic evolution. 

We thank M. Ploskon and A. Poskanzer for discussions 
on the experimental aspects of the proposed observable. 
This work is supported by the Director, Office of Energy 
Research, Office of High Energy and Nuclear Physics, 
Divisions of Nuclear Physics, of the U.S. Department of 
Energy under Contract No. DE-AC02-05CH11231. 



a smaller jet-i^ seems a plausible alternative to viscous 
corrections for a reduced v-i at low to intermediate p t . 

Given the above parametrization we can calculate the 
proposed correlation Cfb using (|10[) . The results are 
shown in Fig(3j While all three cases produce similar 
t>2, they are readily distinguishable by Cfb in the low 
to intermediate p t region. The non-monotonic structure 
seen for scenarios (a) and (b) is due to interplay between 
a non-vanishing jet contribution and rapidly rising hy- 
dro flow contribution to i>2 at low p t , which is absent in 
case (c). In addition, the deviation Cfb from unity at 
low pt provides a measure for the non-collective contri- 
bution to vi in this region and their growth with increas- 
ing pt . While the curves in FigJ3] depend on the specific 
parametrization, they demonstrate the sensitivity of Cfb 
to the non-collectivity which may hide in v-2.. 

In summary, we have proposed to use backward- 
forward correlations of the elliptic anisotropy to distin- 
guish the collective (flow) from non-collective (jet) contri- 
butions to i>2. Using a two-component parametrization, 
we have studied the sensitivity of this observable to the 
degree of non-collectivity. We have further demonstrated 
that this observable is capable of exposing the transition 
from collective flow to jets with increasing p t . 

Concerning an actual measurement, one issue is that 
even though the anisotropy of the hydro flow field is per- 
fectly aligned in the F/B bins, the actual v-i of produced 
particles in the two bins may deviate (both in magnitude 
and in orientation) from the sup posed "U2" due to sta- 
tistical fluctuation (see e.g. [20(] ) . This will reduce the 
correlation from unity even for purely hydro flow, with 
the effect scaling as 1/s/Nun- The other issue relates to 
the low yield in each event from intermediate to high pt 
which may cause the correlation to vanish trivially. These 
problems may be partially cured by selecting proper p t 
bin size and/or event trigger. Also those effects are small 
at low p t region, where the particle abundance is large. 

To conclude, a measurement of Cfb will reveal the 
magnitude of the non-collective contribution to v% and 
the transition pattern of vi between collective and non- 
collective behavior. The centrality dependence of this 
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